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Abstract 

• The transfer matrix of the XXZ open spin-^ chain with general integrable boundary conditions and generic 

anisotropy parameter (q is not a root of unity and \q\ — 1) is diagonalized using the representation theory of 
the q— Onsager algebra. Similarly to the Ising and superintegrable chiral Potts models, the complete spectrum 
I , is expressed in terms of the roots of a characteristic polynomial of degree d = 2 N . The complete family of 
eigenstates are derived in terms of rational functions defined on a discrete support which satisfy a system of 
coupled recurrence relations. In the special case of linear relations between left and right boundary parameters 
for which Bethe-type solutions are known to exist, our analysis provides an alternative derivation of the results 
of Nepomechie et al. and Cao et ai. In the latter case the complete family of eigenvalues and eigenstates splits 
in two sets, each associated with a characteristic polynomial of degree d < 2 N . Numerical checks performed for 
O i| small values of N support the analysis. 

£h ; PACS: 02.30.1k; 11.30.-j; 02.20.Uw; 03.65.Fd 
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^ . 1 Introduction 

Since the pioneering work by Bethe on the Heisenberg isotropic chain [I], various methods and techniques have 
been proposed to solve quantum integrable models on the lattice. Notably, all of them take their roots in 
Bethe's article or Onsager's articles on the planar Ising model [U[3]. By now standard approaches, the so-called 
algebraic and functional Bethe ansatz have been successfully applied to a large number of quantum integrable 
models defined through a transfer matrix or Hamiltonian, for instance quantum spin chains or vertex models with 
periodic boundary conditions. In both approaches, the transfer matrix eigenvalues are written in terms of solutions 
of Bethe equations whereas the construction of the so-called Bethe eigenstates requires either the existence of a 
simple reference state (also called pseudo-vacuum state) on which the creation operator acts iteratively, or the 
construction of modified Baxter— Q operators. In any case, the proof of completeness of the eigenstates is not 
straightforward and makes necessary numerical computations. 

The homogeneous XXZ open spin-^ chain with the most general integrable boundary conditions is among 
the models for which these standard techniques cannot be applied, except for some special family of boundary 
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conditions or special values of the anisotropy parameter. Its Hamiltonian reads 
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where (71,2,3 an d cr± = (oi ± 102)/^ are usual Pauli matrices and A = (q 1 ^ 2 + q~ 1 / 2 )/2 denotes the anisotropy 
parameter. In the present analysis, we will restrict our attention to the (massless) regime — 1 < A < 1 i.e. 
q = exp(</>) with <fi purely imaginary. To exhibit the six independent boundary parameters, below we will sometimes 
use the following parametrization with 9, 9 £ R and a, a e C: 



For the special case of diagonal boundary conditions considered in [H [5] , or for special linear relations between 
the left and right boundary parameters [HIE] (see also [8]) this model has been studied within the algebraic Bethc 
ansatz or functional Bcthe ansatz frameworks. More recently, a Bethe-type solution has also been obtained for 
arbitrary (generic) boundary parameters provided the anisotropy parameter takes special values (q is a root of 
unity) 9 . However, for generic integrable boundary conditions and anisotropy parameter the exact spectrum of 
(Q]) has remained an outstanding problem. Either due to the lack of a pseudo-vacuum state, or due to the lack of 
functional relations of finite order for the transfer matrix, both algebraic and functional Bethe ansatz techniques 
have failed in this regime of parameters. 

In this paper, a solution to this problem is given: we propose the exact solution of the model |T]) for arbitrary 
parameters and q not a root of unity. The approach here presented, initiated in [10] . partly takes its roots in 
Onsager's original work [2J. Mainly based on the representation theory of the infinite dimensional q— Onsager 
algebra (also called tridiagonal algebra) which underlies the integrability of the model [TT], our results for the 
spectrum and eigenstates essentially differ in form from the ones usually obtained using Bcthe ansatz techniques. 
As expected from the method used, they have closed analogy with the ones for the planar Ising and superintegrable 
chiral Potts models [T2l H5] . 

In the next Section, we recall how the transfer matrix of the XXZ open spin chain |T]) can be simply written 
in terms of the generators of the q— deformed analogue of the Onsager algebra and its Abelian subalgebra. Then, 
these generators are shown to possess a block tridiagonal structure in the vector and functional basis proposed 
in [14] for which all coefficients are derived recursively (see Appendix A). In Section 3, the spectral problem 
for the transfer matrix is then considered in details: the eigenstates - written in the vector basis - are derived 
as linear combinations of rational functions of one discrete argument which takes its value within the roots of 
a characteristic polynomial of degree d — 2 N . Then, we consider the special case of linear relations among 
the boundary parameters and discuss the correspondence with each of the two complementary sets of Bethc- 
type solutions obtained in [6] [2 [8] which, all together, cover the full spectrum. For general integrable boundary 
conditions, all the eigenvalues of the transfer matrix are derived from second-order discrete q— difference equations. 
Concluding remarks follow in the last Section. In Appendix B, numerical results for the case of N = 3 sites are 
reported as an example. 

2 The transfer matrix in q— deformed Onsager's basis 
2.1 An alternative description 

In Sklyanin's approach of boundary quantum integrable models, the transfer matrix of the XXZ open spin chain 
is constructed starting from a R— matrix - which statisfies the Yang-Baxter equation - and two K±— matrices - 
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which satisfy the reflection and dual reflection equations, respectively - associated with the two boundaries. In 
particular, from the most general c— number K— matrices derived in [15] one obtains the Hamiltonian (|TJ) among 
the local conserved quantities. However, as pointed out in [10] the transfer matrix can be alternatively written in 
terms of the generators of an Abelian subalgebra of the g— deformed analogue of the Onsager algebra (see details 
below) with no necessary reference to Yang-Baxter-type of algebraic structures. In this point of view, it reads 



N—i 



txxz (u) = J2 Ki+i («) 4+ } i + (u) X {N) (3) 
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((g 1 / 2 ^ 2 + g" 1 /2 M -2) K + ( g l/2 + ? -l/2j K .J e (/ 

+ (g + g" 1 - u 2 - U~ 2 ) N 

({q l / 2 U 2 + g -l/2 u -2 )/t * + ((? l/2 + 9 -l/2 )/t -) e (iV) 

+ (g + g _1 - u 2 - u -2 )^ 

is a family of rational functions in the spectral parameter u. The explicit expressions for the coefficients 

and luq N \ ejf^ essentially depend on the choice of quantum space representations at each site (two-dimensional 

for (TTJ)) of the spin chain. Calculated recursively in [TT] (see also [H]), for the model flj they read 

N 

CW = (-!)*- V/ 2+g -V2 r+ i E «fe-«far— i W 

fel<...<fejV_ n _!=l 

where we have defined 

+ , ^"V 72 - <T 1/2 ) 2 2(g + g^) 

1 (gV2 + g-i/2) + j^fc-feVa+g-i/a) ' afe ( g i/a + ? -i/2j 101 - s .-.^» 

and for arbitrary values of N 

e { ± N) = (g + g- 1 )4 Ar - 1) -2 £ r i ), 

W , lV v fc + fc- / 2(g + g- 1 ) ^-i/ 2(g + g- 1 ) , efeffl (g 1 / 2 - g^ 2 ) 2 X 



1 ' q 1 / 2 -q- 1 / 2 \(q 1 / 2 + q- 1 / 2 )) [(q 1 / 2 + q~ 1 ' 2 ) k+k^q 1 / 2 + q~ 1 / 2 ) )' W 



On the other hand, the algebraic part of for the model (UJ is formed by a family of mutually commuting 
operators acting on the tensor product of N two-dimensional quantum spaces. Generalizing the Onsager's hier- 
archy, they form an Abelian subalgebra [TT] . Indeed, similarly to the undeformed case each operator I^i+i can 
be written in terms of the generators of the g— deformed analogue of the Onsager algebra introduced and studied 
in [H>]. Using the basis of AN operators Wif \ Wi+i , Gi+i . Qi+l with i = 0, 1, JV - 1 introduced in [TT], they 
take the simple formd 

xffr = ^ + + + £C • (6) 

3 Note that the remaining boundary parameters are implicitely contained in the realizations of VVj^> 0jVj? i ■ 
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The fact that all operators 4+1 > ' = 0, 1, JV — 1 are in involution is actually a direct consequence of the 
defining relations (see details below) of the q— deformed analogue of the Onsager algebra [JB]. In view of the 
expansion ^ it follows 

[42i>4+\] = [Hxxz,liV+i] = fOT all fc, i 6 0, .., N - 1 (7) 

which characterize all Abelian symmetries of the XXZ open spin chain JT]) [TJJ . More generally and independently 
of the explicit form of the rational functions J-2Z+1 (u), from ([7]) and ^ one easily recovers the well-known 
integrability condition for (fl]): 

[txxz(u),t X xz(v)] =0 . (8) 

According to ([3]), the spectral problem i.e. finding the eigenstates and eigenvalues of the transfer matrix 
(as well as any local or nonlocal conserved quantities) reduces to the one for the family of mutually commuting 
operators ([BJ. In view of ([BJ, this problem clearly needs a detailed analysis of the finite dimensional representations 
of the q— deformed analogue of the Onsager algebra which is considered below. 

2.2 The q— deformed analogue of the Onsager algebra and tensor product represen- 
tations 

In [TBJ, based on the structure of the Sklyanin operator a q— deformed analogue of the Onsager algebra has been 
introduced and studied in details. With elements W_fc, Wfc + i, Gfc+i, Gk+i, it is a non- Abelian infinite dimensional 
algebra with defining relations 

[Wo,W fc+1 ] = [W_ fe ,Wx] = —L-—^ - G fc+1 ) , 

[W ,G fc+1 ] g - [G fc+1 ,W ] 9 = pW.fe.i-pWfe+i , 
[Gfc+i, Wi] q = [Wi, Gfc+x] q = P \N k+2 - p\N_ k , 
[W 0) W_ fe ] =0 , [Wi,W fc+1 ]=0, 

[Gfc+i, G w ] = , [G fc+1 , G ;+1 ] = , [G fc+ i, G,+i] + [G fc+1 , G ;+1 ] = (9) 

and fixed scalar p, kA £ N, where the q— commutator [X, Y] q = q x l 2 XY — q~ x l 2 YX has been introduced. In 
particular, assumingQ Gi = [Wi,Wo] g and Gi = JWo,Wi] 9 it is not difficult to exhibit a close relationship 
between the algebra ([9]) and the q— Onsager algebr (also called tridiagonal algebra introduced in [17]). Indeed, 
replacing these elements in the fourth relations of ((9]) for k — 1 and using the second and third relations one 
immediately recognizes the so-called q— deformed Dolan-Grady relations discussed in [17] (see also [10]V 

For the XXZ open spin— i chain with generic boundary conditions ([T]), the algebra ^ is a hidden dynamical 

(nonlocal) symmetry of the model: It has been shown that the operators W_; , , G\+1 , Gi+t i n ® induce a 
representation of © provided one identifies [TT] 

p = (q 1/2 + q- 1/2 ?k+k_ . (10) 

By comparison with the case of the Onsager algebra, for finite dimensional representations linear relations among 
the realizations of the elements generating © are well-expected. Indeed, for the finite dimensional representations 
associated with the model (JTJ we have found [TT] 

/ 1/2 —1 /2\ N 



4 It is the case for the corresponding operators of the XXZ open spin chain {T}. 

5 However, an isomorphism between the q— deformed analogue of the Onsager algebra J9} and the q— Onsager/Tridiagonal algebra 



has not been exhibited yet. 
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with nm> ©■ 

Let us finally describe explicitely the structure of the operators , VVj , ^£ , acting on the tensor 

product representation V : (C 2 )® JV . From the results of pj] which were based on the representation theory of the 
U q i/a(sl2) algebra and its connection with the representation theory of (J9j> , one has 



(qV 



(g 1/2 -g- 1/2 ) „ ^^-D^, „ ^^-ia^ %+o 



wffl = (g i/ a+ V^) J8 ^" 1)+ ^ 3/2 ^«r 13 -^m 1) ) 



fc + fc_(q 1 /2 + g -l/2)2 V ^ r 1 1 J (gl/2 + q -l/2)2 < 



^+1 - IM+l ^ J + ^1/2 +g -l/2)2^ 

*m - <» [yi+i y t ) + ^1/2 + ^-1/2)2^ 

+ M<7 - ff-V + 8 + wfff 1 ' - - W^) • (12) 

For the special case I — note that we i dentifyl 



w Wi -n I -0 flWl -flWl _ k+k-jq 1 ' 2 + q 1 / 2 ) 2 (N) 
W l |j=o=0, W_ /+1 |; =0 =0, C/j |i=0-k; U=0 = ( g l/2 _ g -l/2) ^ 



(13) 



and the identity operator M^ N ' = K ® ■ ■ ■ ® I has been introduced. In addition, one has the "initial" c— number 
conditions 

W^=ef, HfW 0) and = <?f = ef eL°V /2 - q~^) . (14) 

Together with (fTTj) , the expressions above determine the action of the nonlocal operators Z^l+i §i ven by © on 
V. Being of a rather complicated form, below we use different basis introduced in [14 for V on which all nonlocal 
operators act in a rather simple manner. 

2.3 Block tridiagonal structure of the q— Onsager generators and their descendents 

As shown in [14], the family of nonlocal operators Wq ,w[ enjoy remarkable properties. For the XXZ open 
spin chain with general boundary conditions |T]), using the definitions (|12p it is actually not difficult to show that 

gf) = , g[ N) = H N \m w ] q ■ (is) 



6 The reader must keep in mind that W ; l J |i =0 + Wl™ J |(=o ,Wi7+ili=o + Wft{\t=o for any N 
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Consequently, using ([9]) one finds that VVq > Wj ^ also induce a representation of the q—Onsageri.e. tridiagonal 
algebra [T7] with defining relations (so-called g— deformed Dolan-Grady relations) and (fTT))) as: 

[wruwrjonj,.] = Hwf'wr 1 ]. (i6) 

In other words, the pair of operators Wq , w[ N ^ is said to be a tridiagonal pair which general definition can be 
found in [17] . Briefly, for generic boundary parameters ([2]) and generic values of q there exists a complete basis 
(resp. dual basis) of V in which Wq N ^ (resp. VV[ ) is represented by a diagonal matrix with degeneracies and 
W[ N ^ (resp. Wq"') is represented by a block tridiagonal matrix. Such basis and dual one have been constructed 
in [14] . so that we refer the reader to this work for details and only recall the main results here. 

Direct basis: Let N be fixed and assume the convenient parametrization of the boundary parameters (|2j). 
Define the family of 2 N possible spins configurations (e^ = ±1,62 = ±1— j£jv_i = ±lj £jv — characterized 
by the quantum numbers n, i such that0 

N , . 

n={N-^ef)/2 i.e. < n < N and i e {1, ...,()} . (17) 
In the canonical basis <S)jLi l = ' = }j on which the nonlocal operators (fT2"|) act, let us introduce the 2 W states 

<J = ® (e e ^ ] £^4 i] ^+*>| +)i + |_ ); ) (18) 
i=i 

with (|17p . For generic boundary parameters and q, these states form a completed basis of V on which the 
tridiagonal pair W^ N \ W[ N) acts H2. Using the explicit expressions (|12p it is possible to exhibit the diagonal 
and block tridiagonal structure of the tridiagonal pair W^ N \ W[ N) in this new basis: 

W^VS = AW^g with AW=cosh(a+(iV-2nW2), 

(?+0 ,„ n , (?) (?-0 

2=1 2—1 2—1 

with n = 0,1,... ,N and j 6 {1,...,(^)}. Note that the explicit form of the coefficients 0&£ , ^n[ij] ■> ^n\ij] ' 
calculated in [T3j - are reported in Appendix A. According to (fTS)) . the nonlocal operators Q[ N ^ and Q[ N ^ also 
possess a block tridiagonal structure in the eigenbasis ([18)) of Wo^" 1 . Plugging (fl9|) in (fl5|) , one finds 

(n+l) (?) 

y l y n \j] ~ / , \1 A " Q A n+l) U n\ii\ + / A1 1 J A « U n[y] V„[j] 

i=l i=l 

(?-0 

+ w A 2i y Vi/^j] v n _i[j] ? 

i=l 

i=l i=l 



7 Given n, the ordering of states choosen here will be described below. 

8 For special relations among the boundary parameters, the representation may become indecomposable. We do not consider such 
possibility here. 
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Proceeding similarly for higher values of I — 1, 2, N — 1, it is an exercise to show that all nonlocal operators 
defined by (|12p also enjoy a block tridiagonal structure in the basis (fT5|) . To show it, it is useful to recall that the 
states (TT5)) are ordered such that 

*sa - (e* f(w - i - aww ^i+^+i->^)®^" 1) for 



n 

( e - a - (N+ .- 2nW2+i e l+)N + h)jv) ^)_ ( _ }] for j e{( N n 1 ) + 1, .... Q 



which are used to obtain recursive formulae for the coefficients of the block tridiagonal matrices. Using ([12]) 
together with the linear relations (fTTj) , after some straightforward calculations we finally get: 

(») 

W -l V n[j] - 2^ a n[ij] V> n[i] , 
i=l 

(? + x) . v „ (?) (?-x) 

i=l i=l i=l 

(n l) (n) 
i=l i = l 

T ^ « 9 »- 1 - 1 U n[ij] Vn-l[i] ' 

2 = 1 

(n+l) (n) 
t=l j=l 

i=l 

where all coefficients are reported in the Appendix A. Note that they have a rather simple form regarding to the 
recursion on the index /, a fact which may have been anticipated from the linear relations pip . 

Dual basis: The previous analysis can also be conducted considering the eigenbasis of w{ introduced in [14j 
instead of Hq , leading to analogous results. Due to their close relationship with the expressions derived above, 
we only report the main results here. Recall that the 2 N states defined by 

form an eigenbasis of w[ N ^ [II]. In this basis, the action of all nonlocal operators reads: 

(" + 0. (jvn (") ("0 

M ,(JV) (AO _ /> W i h( N ' 1 ^ W i W 

m—1 m—1 m—1 

( N ) 

W _ ft ' ' ,(N) 



m—1 



U+0 -rjvn (?) _, 



c co w co _ vrrtVa^w „-i/2 \(w)i a (jv, ° ,„w , v ^ (JV) 

^+l^[fc] - LW A S +1~9 A s ) s[™fc]V s+ i[ m ] + 2^ '^[mfej^lm] 

m—1 



(f-x) 

2—,yl A s-\~1 A s )l'8[mk] i P s -l[ 
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(N \ (N\ 

-(Nl) ^ s I ~ ( N I) 

m=l m— 1 

+ E^^AW-^aWj^^ (22) 



with the identifications 

%i[ml\ = ^mfcjU^-a*, 0^-0 for X £ { (X, 6, C, G?} , 
'^[mfe] — 1 i 8[mk]\a^-a* , 4>^-4> ■ 

Functional basis: Our aim is now to describe the functional basis which provides another useful realization of 
the nonlocal operators wi^' , VV^Tj , G\+l , Qi+} and will play an important role in the derivation of the transfer 
matrix spectrum. Using the scalar product (., .) in the canonical basis <S)jLi | ± ) _? , let us consideiQ the family of 
elements in the discrete argument z s — exp(a* + (N — 2s)(p/2): 



fS^) = (^S^S) with " = 0,1,...,*, i €{!,...,( )} (23) 



where 



N N 

fs[k] _ 

1=1 



iv 

EL"i 4^/2+«| + ) , + with s = (jV- £ 4fl)/2 (24) 



for any s = 0, 1, N and fee {(l, ( )}. Given k, s fixed, there are exactly 2 N elements of this form which 
are easily computed by taking the scalar product of (fT5|) with (|24p. 

Having in mind that (Wq N ^)^ = w[ as <ft — —<j) and — a*, one writes the obvious identitities 

VPs[k]> VV V n [{\) - V V 1 ^s[fc]'^n[i]J anCL l KV ¥> s [fc] ) W n [i\ ) ~ VP s [k\ ' ¥V 1 r n [i] ) * ^ 

Acting with the nonlocal operators Wg^, on the basis (|18p . (|24|) one can show that the family of elements 

{'nh'il ( Zs )} satisfy a system of coupled three-term recurrence relations in n and a system of coupled second-order 
discrete q— difference equations in the argument z s of the form (see |14j for details) 

(n+0 (») (n-l) 

^^iyfe]^ 8 ) = X! F i+i[ 4 fc](^) + E fl «^ F 4i]( z ») + E C ^}? F n-i[ifc](^) > ( 26 ) 

-i— 1 i— 1 2—1 

m— 1 m— 1 m— 1 

respectively. Here, note that we have introduced the "dual" spectrum = (z s + z s _1 )/2. 

As explained in [T3], the equations (f2"o| . ([2T]) generalize the ones defining the Askey- Wilson q— orthogonal 
polynomials on a discret support. The main difference here is the presence of degeneracies in the spectrum of 
Wq N \ W[ N ^: for N > 1, the family {F^ fe j(z s )} is associated with rational symmetric functions of the discrete 

argument As (symmetric under the change z s — > z^ 1 ) with s — 0, 1, ...,N. Associated with each equation in 



Z5]), from (f2"o| , (j2"T|) one identifies the action of Wq , W{ on the space of functions {> n y k ]( z s)}, as well as 

the action of Q[ , G[ N ^ ■ In particular, the operator associated with Wq N ^ can be seen as a generalization of the 
second-order q— difference Askey- Wilson operator [Tl] . 



9 The normalization is different from 1141 . 
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Generalizing the calculations to higher values of I = 1, 2, .., N — 1, it is straightforward to show that any of the 
nonlocal operator Wi^, W^fJ, Gffl, G~i+i acts on the functional basis (|2"5|) as a second-order discrete g— difference 
operator in the argument z s . We eventually get 

(f +1 )_ CJV n (") GU) 

m—1 m—1 m—1 

(f) 

\ < P a [k]> VV l+l%\j]) ~ L a s[mt] h n [ 3m ]( z S J ; 

m—1 

(f+l) (s) 
m=l m—1 

+ 2^ A s ~ 9 A s -lj C s[rofc]r n [^](,32 s J J 

m—1 

(s + l) (s) 

l^[fc]'^+l^nb]J - 2^ W A s+1 - <1 A s )V a [mk]r n y m ]{.q Zs) + 2^ < l s[mkf n \jm]{ Z s) 

m=l m—1 

+ E (l 1/2 ^-\ ~ ?- 1/2 A^)5S] Fjg,, (?*.) • (28) 

m—1 

Alternatively, it is worth interesting to notice that we also have 

As a consequence, using (j2"0)l it is alternatively possible to derive a dual realization of the nonlocal operators 

* life] 1 



acting on the functional basis F^,, (z s ). Although we will use it in the next Section, we leave it as an exercise 



for the reader. 

Now, having in hand all useful realizations (|20p . (f2"2")l . (|28p - and its dual version just mentionned - of the 
nonlocal operators Wl^ , i i i we now turn to the block tridiagonal structure of the transfer matrix 
in each basis. 

2.4 Block tridiagonal structure of the transfer matrix 

From ([5]), the action of the nonlocal mutually commuting operators I^f+i m above basis readily follows. For 
generic boundary conditions, none of the entries of the block tridiagonal matrices are vanishing and all states 
{V>^q} are linearly independent. Then, 

(sr+o (?) (?-x) 

X 2/+l^n[j] — 2-^, n[ij] Vn+l[i] + 2_^ n[ij] V n [i] + 2_^ n[ij] "n-l[i] ' V ou ; 
i=l i=l %=1 

&0 (™) C-0 

(<5 W J w ^ (JV) ) - V R W) F (N) (n- l z ) + V A (N - l) F W fz ) + V C {N ' l) F W (az ) (31) 

771=1 771=1 m — 1 

where the coefHcients read 

sJS? = («* + £(<W> - ?- 1/2 + ^(. 1/2 aS - t 1/2 a^)) , 
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and 

4mi] = 4™l]U K . ; ^Z;A^A for Xe{a,b,C,d,h} and X G {A, B,C} . 
Alternatively, from the action of all nonlocal operators on the dual basis (piTj) one immediately gets 

C+0 (™) ("-i) 

4+w$ = y ®\ N tw ( % i + y A N tWf\ + y c\ n ^ {n ], , , (33) 

2/+lr's[feJ / / s[mfc] r s-|-l[m] / / s|mfc| ^s[mj / / s[m/ej ^ s— l[m\ ' V / 

m— 1 m— 1 m— 1 

(»+l) (n) (n-l) 

$n\j]^2l+lVs[k]) = ^n[»7] ^n4-l [»m] ) + A^-i] ^ ^n[»m] ( Zs ) + ^n-l[im] (^*) ' 

z— 1 i— 1 i=l 

where the "dual" coefficients read 



s[mfe] 


■ (■ 


C 'W) 

s[mfc] 


■ (- 


s[mfc] 





7 W) 
s[mfc] ' 



and 



^^yl^^xi^ fOT xe{a,&,c,rf,/i} and X€{A,B,C}. 



[ij] n[ij] Ik*<->k;^C->2^A— >A 

This demonstrates the remarkable block tridiagonal structure of the transfer matrix (J3J) in each basis. In 
particular, from the quantity (^ffi, txxz(u)tp^J) (or alternatively (^S,txxz(w)^Sj)) one extracts a gener- 
ating function for a family of mutually commuting second-order discrete q— difference operators generalizing the 
Askey- Wilson one. 



3 Eigenstates and spectrum of the XXZ open spin chain 

3.1 Generic integrable boundary conditions 

For generic values of the boundary parameters and deformation parameter q, none of the coefficients (|32p are 



vanishing. In other words, there is no nontrivial subspace of V which is left invariant under the action of the 
nonlocal commuting operators l^i+i- Then, any common eigenstate of all nonlocal commuting operators I^i+i 
can be written either in the eigenbasis of Wq N ^ - i.e. the direct basis (fT5f - or in the eigenbasis of wj 7 ^ - i.e. the 
dual basis (|2Tjl . Choosing for instance the most general linear combinations of states (fTS|) any eigenstate of I^l+i 
can be written 



A 



n=0 j=l 



where A^ j is an argument used to classify all posssible eigenstates of the transfer matrix ([3]). In view of Q and 
the block tridiagonal structure of (resp. Wq N ^ in the dual basis) with respect to the eigenbasis of Wq N ^ 

(resp. W[ N) ) for generic boundary parameters, it is actually sufficient^ to consider the action of x[ N ^ on to 
determine all functions (A^ ) as wen as the discrete support which A^ belongs to. Using ([50)1 for I = 0, 
the 2 N functions are found to satisfy the system of coupled three-term recurrence relations on n (see also [T5]): 

(n-l) („ + l) (n) 

y. B„_' 1 y m ]/ 7l _ 1 [ m ](A 1 ^) + ^ ^„+'l[j m ]/i+l[ m ](^l ) + £ (Ai[jm] _ ^1 ^] m )fi[m] (M ) = ^ (^^) 
m— 1 m— 1 m—1 



10 For generic values of the parameters, the spectrum of Xj^' does not have degeneracies. 
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for n = 0, 1, N and j = 1, (*) with B^g^ = c£°j jm] = and we have defined 



Appart from the occurence of degeneracies j, m which introduce an additional complexity, such kind of equations 
often arises in the theory of q— orthogonal polynomials, and can be easily solved along the same lines. Normalizing 
/o[i] (A^) = 1, 2 N — 1 equations of ([57)1 determine uniquely the hierarchy of rational functions in the (yet unre- 
stricted) argument Aj . In addition, it is easy to see that the remaining equation coincides with 7?W (A^') = 
where 

V (N) ( A W) = det[x (N) _ A W^j (38) 

is the characteristic polynomial of l[ written in terms of the coefficients (|32|) . It determines the set of "allowed" 
eigenvalues!"]: 

pW(A« ff »)=0 Af» e {A|y,..,A5}. (39) 

Consequently, the complete set of 2 N eigenstates of the XXZ open spin chain ((T|) with generic boundary 
parameters and q is given by (f3T>)) with (f37|) . (|3"9")l . Note that the completeness of the basis (j3l))) is a natural 
consequence of the dimension and irreducibility property of V for generic boundary parameters. In our approach, 
the counting of eigenstates is encoded in the degree d of the characteristic polynomial (138|) . For generic boundary 
conditions, the irreducibility of the representations (fT5|) . (|36p implies al = 2 N . 

Let us now turn to the calculation of the eigenvalues of all nonlocal commuting operators ^ffi ■ The general 
spectral problem reads: 

j(f ) *W - A (JV) ¥ N) for / - 1 AT — 1 

X 2i+1* A (N) — Ji 2(+1* A (N) 101 t— U,l,...,iV 1. 

Taking for instance^ the scalar product of each side of the equation with and using (|3~Tj) , we end up with 
a set of second-order discrete q— difference equations on s which - independently of each other - can be used to 
determine the eigenvalues. These constraints, satisfied for any s — 0, 1, N and k = 1, ( JV ), read 

m=l w [fe] W m=l *[fe] l s i m=l I s J 

where the 2 N — different basis of functions 

N (?) 

n=0 j=l 

with (|23p have been introduced. For generic parameters, it turns out that the eigenvalues Ajxz(") of the transfer 
matrix ([3]) are rational functions of the spectral parameter u and A^ defined on the discrete support (|39| . They 
read 

N-l 

Axxz{u) = ^2i+i(u) Aj\ + F (u) . (42) 

1=0 



11 The spectral problem for the trivial special case N = 1 of the XXZ open spin chain JTJ has some analogies with the one for 
the Azbel-Hofstadter model. From J37) one immediatly finds (A* 1 ' ) = 1 , fm](Ai^) = -(^ojii] ~ A i 1> )/ C i["ii] and thc 

characteristic polynomial P^^A^ ) = (-^qJii] — A l ) (^[m] — ^-l) — ®o[llpl[ll] wn ^ cn ^ as exactly two distinct roots. 

12 Alternatively, on may simply consider each of the recurrence relations associated with I > 1 generalizing J3TD for a given value of 
n. In this case, the eigenvalues are written in terms of the coefficients i l'{2[ ) and ratios of /^yj (^i^) - 
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The spectrum E of the Hamiltonian fl} immediatly follows: 

_ g-l/2)(gl/2 + g-Vy ^ dF 2t + 1 ( M ) (iV) AFqW, 



fa 1 / 2 - o- 1 / 2 ! 2. 



2(/c + «*)(e^ + eL u; ) 

(g l/2_ g -l/2)2 

2(g 1 /2 + g -i/2) 

3.2 Special integrable boundary conditions and Bethe-type solutions 

For special relations between left and right boundary parameters, invariant nontrivial subspaces denoted W± C V 
with dim(W±) < dim(V) and W± C V dim(W±) < with dim(V) such that 

dim(W+) + dim(W_) = dim(V) = 2 N , 

dim(W+) +dim(W_) = dim(V) = 2 N (44) 

can be constructed. Indeed, due to the block tridiagonal structure (j3"0)) or (|33p of all nonlocal operators with 
respect to the eigenbasis of VVq and w[ N ^ for some vanishing entries in (f3"2")) or ([33| it is possible to construct 
eigenstates of ^i+i as truncated - compared with (|36[) - linear combinations of states (fT8|) or (j2Tj) . Then, as we 
now show given a family of boundary conditions the complete set of eigenstates of the model ([1} split in two 
sectors: 

Boundary conditions associated with W±: Given an integer < p < N fixed, consider the two families of 
states written in the direct basis (fT5|) and dual basis (|2"Tj) . respectively, as: 



N-p-l 



n=0 j=l 1 s=0 fe=l 



According to ([50)1 and (f3"3")> . the states (|4"5|) can span the two invariant subspaces denoted W± provided the 
conditions 

B pS]°= and S JV-i-i[mfc] = fora11 hi,3,m,k 

are satisfied, respectively. Written in terms of the boundary parameters and using the parametrization |(5J) in 
(|32|) . ([35]) , it is easy to see that both conditions reduce to 

a±& = -i(6-6)-(N -2p- 1)0/2 mod(2in) . (46) 

In the first sector - associated with >V+ -, the family of rational functions f^j\ + \^-i) must satisfy the truncated 
(compared to (|3"T|0 system of three-term recurrence relations 

(n-l) (n+l) (?) 

E ®n-l[jm]fn-l[m}{^ ) + E ^+1W^+1h(^ ) + E C^i^m] ~ M ^i m )/l[mf ' (^1 = ( 47 ) 
m—1 m=l m—1 

with n = 0, 1, On the other hand, in the second sector - associated with W_ - one has 

( - ) ( N ) (") 

E ^^(ah + e <5l]/Sri]( A ^) + E (4f» - A^M/SrW 3 ) - o ( 4 8) 

m—1 m—1 m=l 

with s = 0, 1, N — p — 1. Similarly to the case of generic boundary parameters, dim(VV±) — 1 equations of (|4"T)l 
or (|48[) determine the functions (A[ ) , whereas each of the two remaining equations coincides - up to an 

overall factor - with two conditions on some truncated characteristic polynomials 

V^(A^) = <=> A« e {Ag), ...,AW m(w+) } , (49) 

^(AH = A™ E {A^ _ dlm(w _ )+1 , ...,AS»} (50) 
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where 

V (N,±) (A (A0) = det[T^ c - A[ N) J(dim(w ±))] (51) 

and the two truncated block tridiagonal matrices ^j^-unc which entries appear in (|4"?1) and respectively, have 
been introduced. 

Each linearly independent set of eigenstates in (|4"S")) form an eigenbasis for l{^ c of W± which dimensions!*^!, 
respectively, follow from the degeneracies of each eigenspace of W^ N) and W[ N) One has 

dim(W+) = ^( J and dim(W_) = ( J (52) 

which clearly satisfy (|44)) . Then, given the family of boundary conditions (|46[) with < p < AT the complete set 
of eigenstates (±) of the transfer matrix ([3]) is given by (|45[) with (|49f for (+) and (|50| for (— ). 

Boundary conditions associated with W±: Given an integer < p < N fixed, let us consider instead the two 
families of states defined by: 



n=p+l j=l s =N—p k=l 



*iW= E E/S +) (ArX-J and ^7 3 = E E/Sr } (ArVS- (53) 



Following the same procedure as above, it is easy to see that these states can span the invariant subspaces denoted 
W± provided 

C^'J). = and cS r N ' 1 } ,,=0 for all l,i,j,m,k . 

p+l[tj\ N—p[mk\ ' SJ ' ' 

In terms of the boundary parameters, these constraints become: 

a ±& = i(9-6)-(N-2p- 1)0/2 mod(2m) . (54) 

The corresponding families of rational functions entering in (|55|) satisfy the system of truncated three-term 
recurrence relations (|47| with n = p + 1, N and (|48|) with s = N — p, N, respectively. From these latter 
relations, the two characteristic polynomials 

V^ ±] (AH = det^tl - (55) 

are easily exhibited. Here, by analogy with the previous case we have introduced the two truncated block 
tridiagonal matrices I^trunc- For each family of states (±) in (|53|) . it turns out that the eigenvalues belong 
to a discret support defined by: 

?<"<><»>) =0 «. A<»>e(A<»>,...,A« m( _ )} , (56) 

^-'(AS^-O - AS»',(A«_ dim , w _, + ,^A<»>,^ (57) 



In view of ([53]) . note that one has (l44|) with 



dim(W+)= E ( ) and dim(W_)= ^ ( J . (58) 

n=p+l ^ ' s=N-p ^ ' 

Then, given the family of boundary conditions (j54"|) with < p < N the complete set of eigenstates (±) of the 
transfer matrix ([3|) is given by (|53| with (|56f for (+) and (|57|) for (— ). 



13 Similar results have been independently obtained from the representation theory of the two-boundary Temberley-Lieb algebra 

ED- 
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The complete set of eigenvalues of the transfer matrix ([3]) can now be calculated. For any choice of boundary 
conditions (I46[) or ([51)) , as a consequence of the eigenspaces decomposition the eigenvalues of all nonlocal operators 
^2/+i s P u t m two sets. Introducing the truncated functions 

= and *W->(») = , (59) 

by analogy with the case of generic boundary conditions using pip and (|34|) we end up with two independent 
sets of second-order discrete q— difference equations on s and n which determine the eigenvalues for any s,n, k,j. 
We obtain 

AN,+) _ -ST kW) _M y S + 1 ) \ - 7(jV,i) [to] ^) M I s ^ , fim 

A 2(+l - 2v °«[mfc] v,v(W,+)/x + ^[mfe] (AT, + ), + ^[mfe] ^(JV.+J/x ' ^ 

m=l *[fc] l s J m=l *[fc] \ s ) m=l w [fc] l s i 

m=l W m=l W m=l \ n > 



All together, using (|42[) they provide the complete set of eigenvalues of the transfer matrix as well as the energy 
levels (|43|) of the model ([1]) for the family of boundary conditions (|46|) , (|54[) . 



It is worth interesting to notice that for the boundary parameter spaces ([54")) , the transfer matrix ([3]) is 
known to admit Bethe-type solutions [Si [3 [8]. Indeed, acting with the gauge-transformed creation/annihilation 
operators B(u), C(u) on two different pseudo-reference states (see [7J[H] for details) the complete set of Bcthc 
eigenstates are formally written 

V(u 1 ,...,u p )=B(u 1 )...B(u p )\tl + ) and 9(u x , ...,u p ) = C(u 1 )...C(u p )\Q-) (62) 

where {uk} denote the Bethe roots. From our analysis, given the family of boundary conditions (|46|1 or (|54| it 
follows that the eigenstates ([521 admit an alternative description of the form ([4"5[) or ([55]) , respectively, as rational 
functions of the discrete variable A^. In particular, setting p = in (|46[) . (f54[) one identifies the two reference 
states 

|fi + >-^$ and |n-)-^^ (63) 

in agreement with [71 [8] . For generic values of p and boundary conditions of the form (|4"6"[) , ([54")) the correspondence 
between (j62[) and ([45]) or ([53]) allows to derive exact relations between Bethe roots and roots of the characteristic 
polynomials V^ N,i '\K^),'P (A[ N ^). These relations will be considered elsewhere. 



4 Concluding remarks 

The approach initiated in [10] and based on the representation theory of the q— Onsager algebra opens the 
possibility of studying massive continuum or lattice quantum integrable models from a point of view which 
possesses some advantages compared to standard - like algebraic or functional Bethe ansatz - approaches: analytic 
solution for generic boundary and anisotropy parameters; no ansatz required; rigorous proof of completeness. 
Applied here to the XXZ open spin chain it has shown its efficiency for generic integrable boundary conditions 
in ([1]) as well as for the special conditions (|46f . ([54[) for which the complete spectrum has been derived. No 
modifications of the formalism between the two regimes of boundary parameters have been necessary, the only 
changes being the structure of the eigenstates and related characteristic polynomial. 

Another nice feature of the present approach relies on the fact that the eigenvalues of the transfer matrix 
for any integrable boundary conditions are expressed in terms of the roots of a characteristic polynomial of 
degree d = 2 N (d < 2 N for ([46]) . ([54]) which gives all subsets of the eigenvalues), similarly to the Ising [2] and 
superintegrable chiral Potts models [T2], [13] . For the regime of parameters for which the Bethe ansatz applies, one 
is reduced to solve a single algebraic equation instead of highly transcendental equations (Bethe equations) in the 
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Bethe ansatz approach. Numerically, a detailed comparison between both approaches is an interesting problem, 
which however goes beyond the scope of this paper. In particular, in the algebraic sectors (|46[) . (|54[) associated 
with Bethe-type solutions, it is worth noticing that such a correspondence between these roots and solutions of 
Bethe equations should not be surprising: it is a well-known fact in the context of q— orthogonal polynomials 
theory. 

Although wc didn't consider such possibilities here, for q a root of unity, or dynamical integrable boundary 
conditions^, or the presence of inhomogeneities and higher spin representations located at each sites, it is clear 
that our approach can be applied analogously: in this case cyclic or higher tensor product finite dimensional 
representations of the tridiagonal algebra need to be considered in details (see e.g. [TT]). In this direction, infinite 
dimensional representations of the tridiagonal algebra also clearly need further investigations: they should play 
a crucial role in the analysis of the thermodynamical properties of the spin chain (N — > oo) for generic boundary 
and anisotropy parameters, still unknown. 

To conclude, using the complete basis of eigenstates (j56")l . (|45[) or (|S"5)) the problem of calculating correlation 
functions in the XXZ open spin chain with integrable non-diagonal boundary conditions ([T]) can be handled. This 
will be the subject of a separate work. 
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Appendix A: Coefficients of the block tridiagonal matrices 

All coefficients below have been calculated using the action of the nonlocal operators as defined in (TT2|) on the 
family of states (IS]). For 1 = 0,1, ...,N - 1 we find: 

Coefficients 

For z^eU,...,^- 1 )}: 

t N ,l) _ COsh0 (N,l-1) COsh0 AN -1,1-1) ^ 1 (N-ld-1) 

~ cosh 2 (</>/2) n[ijI cosh(0/2) cosh(0/2) n M 

sinh(q + {N - 2 - 2n)g/2) f (n-u) _ i^-u-i) 
sinh(a + (N - 1 - 2n)<£/2) \ n M U ™M 



(N-l,l-l) 



H 

2 cosh 2 0/2) sinh(a + (N - 1 - 2n)0/2) " M 

For i€{(Z-l)+l,...,(Z)}, , |1 C ')[: 

m) = cosh0 (N,i-i) 2e a +( jV - 2 ")^/ 2 sinh(0/2) cosh 2 (a + (JV - 2n)</>/2) ( N -i,l-x) 
n[ij] cosh 2 (</)/2) n[ij] cosh(0/2) sinh(a + (N + 1 - 2n)0/2) »[<-(?- x ) o] 

, e °+( jv - 2 ")^ 2 sinh(0) c( ^_ 1)0 



sinh(a + (N + 1 - 2n)0/2) i] ' 

For j6{n+l,.,(?)}: 



14 For integrable dynamical boundary conditions, the corresponding open spin chain also enjoys a hidden q— Onsager dynamical 
symmetry |20| . 
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W ) = cosh0 2e-"^( Af - 2 ")^ 2 sinh(0/2) cosh 2 (a + (TV — 2n)0/2) , (at-u-i) 

cosh 2 (0/2) cosh(0/2)sinh(a+(/V-l-2n)0/2) ""^^'-l" ')> 

e -«-(JV-2„)0/2 sinh( -^ ,(JV-1,0 

~sinh(a+(A^-l-2n)0/2) "-![< ' 
For i,j € {ft- 1 ) +1,. ..,$)}: 

(N,l) _ COSh0 (N,l-1) , COsh0 7(JV-l,i-l) 1 „(JV-l,i-l) 

"Wl " cosh 2 (0/2) "Wl cosh(0/2) - 1 ) M» -1 )] cosh(0/2) n-ip-^" 1 ) j-^ 1 )] 

sinh(a + (iV + 2 - 2n)0/2) / (jv-i.i) rfJV-i,l-i) 
+ sinh( a +(/V + l-2n)0/2) ^n-ili-^" 1 ) i-^" 1 )] ""-iM?" 1 ) M?" 1 )] 

1 ,(JV-1,J-1) 



2cosli (0/2) sinh(a + (TV + 1 - 2n)0/2) 



(64) 



Coefficients 0„[y] : 

For iea.-X+i 1 )}' ./• !1 0* ')[: 

cosh0 7(iv,;-i) , /2 ( h {N-i,i) ,(jv-i,;-i)\ 
n[y] " cosh 2 (0/2) + 6 I ~ J • 

For iEO+1 

.(AM) _ cosh0 Ba+(JV _ 1 _ 2wW2 sinh(0/2) f„(JV-i,l) ^-^''-^ ^ 

nfel " cosh 2 (0/2) + sinh(a + (TV - 1 - 2n)0/2) i] ^M^ 1 ) d J 

g+(JV-l-2nW2 1 COsh(o; + (AT — 1 — 2n)(j>/2) i(N-l,l-l) 

2cosh 2 (0/2) sinh(a + (AT — 1 — 2n)0/2) "M"^ 1 ) d ' 

For i€{l P}, iem+l-C)}: 

7 (AM) _ cosh0 tW-i) 
" W1 " cosh 2 (0/2) ' 

For i e + 1, .... i e {tf" 1 ) + 1, .... (?)}: 

7 (AM) cosh0 i(na-i) _, /2 smh{a + {N + l-2n)4>/2) f i(N-i,i) ,(jv-i,j-i) 
^ = C osh 2 (0/2) + 6 sinh( a + (TV - 1 - 2n)0/2) {^H^) ^Mik 1 ) 

Coefficients 

For iell,...,^ 1 )}, 

(jv,o = cosh0 (jy,;-i) 0/2 sinh(g + (-/V - 1 - 2n)0/2) / (jv-i,;) _ (jv-i,i-i)\ 
n[l3] cosh 2 (0/2) sinh(a + (N + 1 - 2n)0/2) V n M n M ) 

For ie{(^_- 1 1 )+l,...,(^ 1 )} ) ie{l,...,(^- 1 )}: 
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(N,l) = COSh0 (N,l-1) 

n[ij] cosh 2 (0/2) 

For i€ {!,..., (£?)}, je^+l, ...,(*)}: 

(JV ,o _cosh0_ _ - a -(JV+i-2n)^/2 sinh((/)/2) ( n (N-i,l) _rf N ~ U ~ 1] 

cosh 2 (0/2) ' sinh(a + (iV + 1 - 2n)0/2) ^ n-ifo--^" 1 )] "^M* 

a -(JV + l-2n)^/2 1 cosh(Q: + (JV + l-2n)0/2) »(Ar-i,t-i) 

2 cosh 2 (0/2) sinh(a + (N + 1 - 2n)0/2) "-^M" )1 ' 

For ^{Cr/l+l,..,^)}, J€{n + 1 (?)}= 

(jv,o _ cosh0 (jy,;-i) 0/ 2 ( r (N-i,i) _ (;v-i,;-i) 

" cosh 2 (0/2) + r-iI-Cr/) M» J -(^ r - 1 )] 



Coefficients <* n ryi : 

For i,je{l (H}: 

iW) C osh0 iW-i) sinh 2 (0/2) (JV _i ; ;_i) _ 1 j N -i,i-i) 
nM ~ cosh 2 (0/2) cosh(0/2) "»[«] cosh(0/2) " M 

+ cosh ( 0/2)Ci u) - sinh(Q+(iV ' 1 :f n)0/2) Ci u '- ) 

M 2 cosh 2 (0/2) " M 

For ietf?" 1 ) +!,...,(?)}, J6{l,...,n}: 

jW) COsh0 W-l) a+ (AT^2n)^/ 2 S inh ('/') sinh ( a + ~ 1 ~ 2n)0/2) 

n[ij] ~ cosh 2 (0/2) " 6 cosh 2 (0/2) ^-(n- 1 )^]' 

For <6{l,...,(y- 1 )},i€{(y- 1 ) + l,...,(jr)}: 

,(JV,0 cosh0 JNA-i) _ a _ (A r_ 2n) /2 sinh(0) sinh(a + (N + 1 - 2n)0/2) ,(Af-i,t-i) 
^ = cosh 2 (0/2)^ +e c^h^) J-CiT- 1 )] • 

For i,je{(^)+l,...,(^)}: 

i(JV,o _ cosh0 tW-i) sinh 2 (0/2) (JV -i,i-i) 1 *n-i,i-i) 

" cosh 2 (0/2) + cosh(0/2) a n-i[i-(--) - CO sh(0/2) fl "- 1 [-(^ 1 ) M?" 1 )] 

, , ./on^- 1 .') sinh(a + (TV + 1 - 2n)0/2) , (Ar-i,i-i) 
+ cosh(0/2)rf„_ 1[i _ (nN - 1)j _ (ri)] + 2cosh 2 (0/2) ^-ili-^" 1 ) i-^" 1 )] • 



Coefficients »l n [y] : 

For {!,..., (^- 1 )}: 



7 (Jv,;) _ cosh0 h (N,i-i) , (jv-1,0 _ , (jv-i,;-i) 
~~ cosh 2 (0/2) fol fol fol 



sinh0sinh(0/2) 



sinh(a + (TV — 1 - 2n)0/2) 



(iv-i,;) „(;v-i,;-i) j(JV-i.O 

tt n[ij] - tt n[ij] + U n[ij] ~ J • 
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For ie + je{i,...,n^ 

7 (JV,0 _ cosh</> i(N.i-i) 
n[ii] "cosh 2 (0/2)^ fol 



V yV/ , smh(a+(N+l-2n)(f>/2) \ 



_ r (N-l,l-l) 



For ie{l ) - ) r)} ) J6{n+l 1 - ) ( n W )}: 
7 W) _ cosh0 7 (Jv,;-i) 
3 ' ~cosh 2 (^/2)^ Iyl 



o -c-W-2nU/2 • v , • utiles COsh(« + (iV - 2n)(j)/2) (i (N-l,l) WJV-X.i-l) \ 

-2e (- 2 »)^smh 0^(0/2) ^ - ^ _ i _ ^ (A- lfo --(^)] ~ K^-^ ) 



For iJe^-^+l, ...,(*)}: 

h (N,l) _ COSb.0 h (N,l-l) , (AT-1,0 ,(JV-1,Z-1) 

fol ~ cosh 2 (0/2) fol "'''n-ifo'] 
sinh 0sinh(0/2) 



sinh(a + (TV + 1 - 2n)0/2) 

( n {N-l.l) (A r_i,;_i) J(JV-1,0 ,(JV-1,Z-1) \ 



The "initial" conditions for any N are given by: 



( jv,-i) _ t(*,-i) _ (jv,-i) _ jW-i) _ 



= tf /2 - ^ 1/2 )^ } Cf . - -(^ /2 - <T 1/a ) eosh 2 (0/2) 5 , 

and 



(i,o) cosha* sinh(a — 0/2) — sinh(0/2) n,o) cosha* sinh(a + 0/2) + sinh(0/2) 

a °[ n l shm^ ' = shrn^ ' 

7(1,0) e Q (cosha + cosha*) sinh(0/2) ( X 0) e _a (cosh a + cosha*) sinh(0/2) 

°o[ni = ' c i : " 



sinh a sinh 



a 



Appendix B: Numerical results for iV = 3 sites 

For small values of N, arbitrary or special values of the boundary parameters and q generic in JT]) we have 
checked numerically (using MAPLE) the block tridiagonal structure (|30|) , (pTT| , ([55]) , of all nonlocal commuting 
operators. We have also checked that the eigenvalues obtained by direct diagonalization of ^ with (fT2"|) coincide 
with the ones derived from (jit))) . The same computations have been done for the Hamiltonian: it shows perfect 
agreement between the eigenvalues derived from (|T|) or (|43p . Being technically reduced to the calculation of the 
roots of a characteristic polynomial of degree d = 2 N (d = dim(W±) < 2 N or d = dim(W±) < 2 N in the truncated 
(algebraic) sector associated with Bethe-type solutions) and eigenvectors of the block tridiagonal matrix l[ N ^ (or 
truncated ones l[ N ^\trunc, l\_ \trunc), the method can be easily extended to higher values of N. 
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Below we report the results for TV = 3 as an example. On the right hand side of the tables, for generic boundary 
conditions the energy levels have been computed replacing (HP)) into (|4"3")l . For special boundary conditions, the 
eigenvalues corresponding to the Bethe eigenstates written in the form ^^n) ^ are underlined. They have been 

computed by plugging (fBTj) into (|4"3")l . The remaining eigenvalues associated with the Bethe eigenstates 

have been similarly computed using the complementary subset <[60j) . Note that for p = (N — l)/2 = 1, one of the 
solutions to (|5"0)) and ([57)) yields the ground state energy as expected from the results in [3] . 

• Arbitrary boundary parameters: For a — 1.3 + 0.7i, tf> — 1.2i, 9 = 0.1, 9 — 0.21, k = 1.2 + O.lli: 



Eigenvalues 



Hxxz © 


Hxxz (O 


2.970995654 
2.198656210 
1.602395044 
0.946000118 
-0.019184890 
-0.055723904 
-3.180523377 
-4.462614856 


2.970995806 
2.198656341 
1.602395091 
0.946000104 
-0.019184648 
-0.055723820 
-3.180523355 
-4.462614866 



• Special boundary parameters: For a ± a = —i(9 — 9) — (N — 2p — 1)^/2 mod(2iir) , 

a = 1.3 + 0.7i, <j> = l.2i, 9 = 0.1, 9 = 0.21: 



Eigenvalues for p = 



Hxxz © 


Hxxz & 


3.618207615 
2.779649839 
1.727013687 
0.001649161 
-0.026778318 
-0.198835568 
-2.826978951 
-5.073927468 


3.618207522 
2.779650028 
1.727013757 
0.001648996 
-0.026778318 
-0.198835520 
-2.826979214 
-5.073927520 


Eigenvalues for p = 1 


Hxxz © 


Hxxz & 


2.493378633 
2.412226563 
1.848074870 
1.636051825 
-0.270698292 
-0.304708338 
-3.835675623 
-3.978649644 


2.493378482 
2.412226647 
1.848075023 
1.636051957 
-0.270698432 
-0.304708400 
-3.835675507 
-3.978649645 
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Eigenvalues for p = 2 



Hxxz © 


Hxxz ([43]) 


2.423171820 
2.358362091 
1.963670555 
1.210097504 
-0.100200908 
-0.193998388 
-3.440275504 
-4.220827174 


2.423171828 
2.358362193 
1.963670532 
1.210097648 
-0.100200806 
-0.193998355 
-3.440275543 
-4.220827207 



• Special boundary parameters: For a ± a — i(9 — 9) — (N — 2p — 1)^/2 mod(2iir) , 

a = 1.3 + 0.7i, <p = l.2i, 9 = 0.1, 9 = 0.21: 



Eigenvalues for p = 



Hxxz © 


Hxxz 63} 


4.906122870 
3.889167375 
2.586895038 
0.451087518 
-1.198292710 
-1.359411089 
-3.053875109 
-6.221693900 


4.906122692 
3.889167292 
2.586895406 
0.451087370 
-1.198292642 
-1.359411137 
-3.053875627 
-6.221693900 


Eigenvalues for p = 1 


Hxxz © 


Hxxz (E3|) 


2.388082920 
2.313332883 
1.826977133 
1.655209880 
-0.212564780 
-0.234717115 
-3.807430794 
-3.928890122 


2.388082936 
2.313332774 
1.826977088 
1.655209967 
-0.212564681 
-0.234717212 
-3.807430776 
-3.928890149 



Eigenvalues for p = 2 



Hxxz © 


Hxxz ® 


2.520902297 
2.455383157 
1.907373082 
1.121112406 
-0.101546424 
-0.226017760 
-3.360315576 
-4.316891175 


2.520902258 
2.455383158 
1.907373054 
1.121112526 
-0.101546492 
-0.226017706 
-3.360315615 
-4.316891206 
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